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We derive a factorization formula for the double Drell-Yan cross section in terms of double parton 
distribution functions (dPDFs). Diparton spin and color correlations and parton-exchange interfer- 
ence terms contribute, even for unpolarized beams. Soft radiation effects are nontrivial for the color 
correlation and interference contributions, and are described by non-perturbative soft functions. We 
provide a field-theoretic definition of the quark dPDFs and study some of their basic properties, in- 
cluding discrete symmetries and their interpretation in a non-relativistic quark model. We calculate 
the renormalization group evolution of the quark dPDFs and of the soft functions. The evolution 
receives contributions from both ultraviolet and rapidity divergences. We find that color correlation 
and interference effects are Sudakov suppressed, greatly reducing the number of dPDFs needed to 
describe double parton scattering at high energy experiments. 
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production, P1+P2 —> In some hadronic collisions, 

two partons in one hadron can have simultaneous hard in- 
teractions with two partons from the other hadron. This 
process is called double parton scattering (DPS)qA rep- 
resentation of DPS in space-time is shown in Fig!]l] DPS 
was first considered in Ref. pj[4] and was subsequently 
studied in the context of jet production [SH?], double 
Drell-Yan [8] [9] and W+jets [TU]. Two examples we will 
study are double Drell-Yan, P1+P2 —> t[t~{^~2^2 ( wnere 
t\ and £2 could be the same flavor), and same sign W 
pair production, pi + P2 — > 

DPS is higher twist, i.e. it is suppressed by order 
Aq CD /Q 2 compared to single parton scattering (SPS), 
where Q is the scale of the short-distance interactions. 
Heuristically, this arises because the two partons which 
collide in the second hard interaction have to be within a 
transverse area of order 1/Q 2 of each other, whereas they 
could each be anywhere within the incoming hadrons (of 
transverse area ~ 1/Aq CD ). 

Experimentally, DPS has been studied in four-jet 
events by the AFS collaboration [11] at y/s = 63 GeV, 
the UA2 collaboration [12] at y/s = 630 GeV and the 
Tevatron [T3] . At the Tevatron, DPS has also been stud- 
ied in 7 + 3 jet events [14] [15] and there is an analysis 
using early LHC data for DPS in W + 2 jets [16]. In 
these experiments DPS is quantified using an effective 
cross section cr e ff, denned in Eq. The measured val- 
ues of a e ff range from 5 to 15 mb. 

DPS is an important background for light Higgs 
searches in the channel pp — » WH — >• tvbb [17-20]. A 
clean channel for studying DPS at the LHC is provided 



I. INTRODUCTION 

In high-energy hadronic collisions, one parton from 
each hadron can collide via a hard interaction to produce 
a final state with a large invariant mass. The classic ex- 
ample of such single parton scattering (SPS) is Drell-Yan 



1 We will denote the beams by their momenta pi,2, without spec- 
ifying the hadron. The most common cases are pp collisions 
(Tevatron) or pp collisions (LHC). 

2 At the Tevatron and LHC one needs to be careful to distinguish 
double parton scattering from pile-up, i.e. two single parton scat- 
terings involving different pairs of hadrons during the same bunch 
crossing, since this would produce a similar signal. Separating 
the two relies on identifying the vertex of the hard collision. 
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FIG. 1. Double part on scattering in space-time. The two 
hadrons have transverse size of order 1/Aqcd and longitudi- 
nal size of order 1/(7Aqcd), where 7 ^> 1 is the boost. The 
two hard interactions are shown by black dots. They have 
longitudinal and time separation 1/(7Aqcd), and transverse 
separation 1/Aqcd- 



by same-sign lepton searches, i.e. via pi + p 2 — >• 
with W + H> £+1/, since SPS is suppressed [2TU24]. The 
maximum incoming partonic charge in SPS is +1 from 
ud, so the conservation of electric charge requires the 
presence of at least two additional jets in the final state 
via ud —> W+W+ud, leading to a large suppression of 
SPS by [a/(47r)] 2 [a s /(47r)] 2 relative to single Drell-Yan. 
There is no corresponding suppression of DPS, which can 
proceed via uudd — » After typical cuts the DPS 

cross section is of fb order [23, 24 , making it more of a 
long-term goal at the LHC. 

In the original DPS formalism, the cross section is de- 
scribed as [5] 




FIG. 2. Double PDFs in momentum space (forward diagram). 

The dPDF in momentum space is defined b}Q 

F^(xi,x 2 ,r_L,/i) = J d 2 z±e~ ir± ' z± F ij (x 1 ,x 2 ,z±,LL) .(2) 

It should be clear from the context whether the last ar- 
gument refers to position or momentum. The dPDF in 
momentum space is shown in Fig. [2] The figure also 
shows that the dPDF is not the squared absolute value 
of an amplitude, since the partons on the two sides of the 
cut have different momenta. 

It is commonly assumed that the dependence on the 
transverse separation is uncorrelated with the momen- 
tum fractions or parton flavors, 

Fij(xi,x 2 ,z±,fi) = Fij(xi,x 2 , fj)F(z±, fj) . (3) 

In addition, a factorized ansatz is made 

F ij {x 1 ,x 2 ,n) = fiixt, /j J )f j (x 2 , fi) 0(l-x 1 -x 2 )(l-x 1 -x 2 ) , 

(4) 

where / denotes a single PDF and the last factors 
smoothly impose the kinematic constraint x\ + x 2 < 1. 
For small momentum fractions 1 — x\ — x 2 ~ 1 these fac- 
tors can be neglected, and the cross section in Eq. ([I]) 
becomes 



i,j,k,l 

Xa i k(x 1 X 3 y/s,fl)&ji(x2X4y/s,fl) . 



■ (5) 



(1) 



Each incoming hadron is described by a double parton 
distribution function (dPDF) F and the short-distance 
processes are described by partonic cross sections <r, in 
analogy with SPS. Fij{x\, x 2 , z±) is the number den- 
sity for simultaneously finding two partons with flavors 
i : j = g, u, u, d, . . . , longitudinal momentum fractions 
xi : x 2 and transverse separation z± inside the hadron. 
Our convention is that in formulae such as Eq. ([!]), the 
first (second) dPDF is for the beam with momentum pi 
(p 2 ). The d-ik(xiXsy/s) is the partonic cross section for 
partons z, k going to the desired final state. S is a sym- 
metry factor that can arise if there are identical particles 
in the final state. 



The effective cross-section is a measure of the area of the 
proton ~ 1/Aq CD , consistent with the fact that DPS is 
higher twist. 

The expression in Eq. is an example of a fac- 
torization formula, where the hadronic cross-section is 
written as the convolution of hard-scattering partonic 
cross-sections which are target-independent, and non- 
perturbative distribution functions which depend on the 
hadronic target. We will see that there are several impor- 
tant modifications to the naive expressions above. In par- 
ticular, the product forms of Eqs. Q and Q are spoiled 



3 We follow the convention of Refs. |25l I26| , where r_L is a trans- 
verse momentum, not a coordinate. 
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FIG. 3. Extraction of the double parton scattering contribu- 
tion to pp — »• W + 2 jets by the Atlas collaboration 16 . The 
extraction is performed by comparing the observed spectrum 
of the normalized total transverse momentum of the two jets 
Afets = |pi_L +P2_l|/(|pi_l| + |P2_l|) to templates for SPS and 
DPS obtained from Monte Carlo programs. (ATLAS Experi- 
ment © 2012 CERN.) 



FIG. 4. Mixing between the gluon single PDF and the qq 
double PDF. 



by QCD radiative corrections. There are also several 
different spin and color structures that enter the factor- 
ization formula, even for unpolarized beams [25H27] . 





In processes such as double Drell-Yan, one can mea- 
sure the invariant mass and transverse momentum of each 
lepton pair. Double Drell-Yan can arise from a double- 
parton process, (qq — >• 7*) + (qq — >• 7*). The transverse 
momentum of each incoming parton is the intrinsic trans- 
verse momentum of a quark in a hadron, and is typically 
of order Aqcd- Double Drell-Yan can also arise from a 
single parton process such as (g — » qq) + (g — » qq) where 
a quark from one gluon and an antiquark from the other 
gluon annihilate into a virtual photon. In this case, the 
transverse momentum of the quarks (and thus the lep- 
ton pair) can be of order the hard scale Q. There is an 
inevitable overlap between the SPS and DPS contribu- 
tions to a physical process such as double Drell-Yan in 
the region where the transverse momentum of the lepton 
pair (not the individual leptons) is small. The total SPS 
contribution to double Drell-Yan is leading twist, but the 
contribution to the small transverse momentum region is 
Aq CD /(5 2 suppressed (because it is a fraction Aq CD /(5 2 
of the total phase space) and is the same order as DPS 
in this region. 

In light of the above discussion, observables that are 
sensitive to the regions of phase space where SPS and 
DPS are of the same order, are used to separate SPS 
and DPS contributions [28]. In a recent analysis by the 
Atlas collaboration [16], DPS was studied in pp — » W + 2 
jets, for a leptonically decaying W. The normalized total 
transverse momentum of the two jets, A™ ts = |pi_L + 
P2_l |/(|pi_l I + |P2_l|) ? was used to separate DPS from 
SPS, and is shown in Fig. [3] 

The dPDF is a new non-pert urbative object. It 
has a renormalization group (RG) evolution similar to 
that for the conventional (single) PDF. The evolution 
of Fij (xi , x 2 , fi) was determined a long time ago in 
Refs. [22l EDI- It has recently been extended to include 
the zi dependence [25j[26]. The RG evolution for the qq 
dPDF is 



a s C F 



Pqi{^)Fiq(z,X 2 ,Y^ll) + P gi ) , Z, r ± , /l) 



/ Xl \ fg(X 1 +X 2 ,ll) \ 

1 \x 1 + X 2 J X\+X 2 J 



(6) 



where P denotes the usual PDF splitting functions. The 
first two terms describe the independent evolution of each 
of the partons with the standard (single) PDF kernel. 
The second line is the contribution of the single gluon 
PDF feeding into the dPDF via g — >• qq, shown in Fig. [4] 
This term has no dependence, and leads to a inde- 
pendent contribution to F(xi, x 2 , rj_, /i), or equivalently 
to a S^ 2 \z±) contribution to F(xi, x 2 , z_l, /i). This 5- 
function in transverse position space leads to a divergence 
in the cross section Eq. ([!]), as discussed in Refs. [25) 126]. 
The resolution of this singularity is related to the issue of 



double counting between SPS and DPS mentioned ear- 
lier, and will be discussed in a subsequent paper [31 . The 
double-counting problem only enters through the mixing 
between single and double PDFs, which we therefore also 
postpone to Ref. [31 . The factorized form in Eq. Q for 
the dPDF is not preserved by the evolution [32] [33] . Sum 
rules for dPDFs were derived and used to find a new 
ansatz that satisfies the evolution equation reasonably 
well [34]. 

In DPS, the two partons extracted out of the proton 
can be correlated in spin and color, as was first dis- 
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FIG. 5. Interference contribution to DPS. Note the orienta- 
tion of the fermion lines. 



to the well-known factorization theorem for single Drell- 
Yan. In contrast to Eq. ([!]), we will include spin and color 
correlations and interference contributions, which natu- 
rally arise. In Sec. |IV| we define the quark dPDF, classify 
its spin and color structures, study its properties under 
discrete symmetries and give an interpretation in the con- 
text of a quark model. The calculation of the RGE for 
the quark dPDFs and the corresponding soft functions 
is given in Sec. [V] It includes a brief introduction to the 
topic of rapidity divergences and the technology of the ra- 
pidity renormalization group with explicit examples. All 
spin and color correlations as well as interference effects 
are considered. We conclude in Sec. ED 



cussed in Ref. [27]. These correlations are not present 
in phenomenological models and were recently revisited 
in Refs. [25] [26] for transverse momentum dependent 
dPDFs. 

In this paper our main focus will be on formulating 
the QCD factorization theorem for double Drell-Yan and 
related processes in terms of dPDFs, and on studying the 
RG evolution of color correlations in quark dPDFs. We 
find that the color-correlated dPDFs contain rapidity di- 
vergences, which are tied to the presence of large rapidity 
logarithms of r 2 ± /Q 2 . We will treat the rapidity diver- 
gences and resum the corresponding series of logarithms 
using the recently introduced rapidity renormalization 
group [35 , 36 . We find that the effects of soft gluon ex- 
change do not cancel for the color-correlated dPDF: the 
color-correlated dPDFs contribute to the cross-section in 
combination with a non-perturbative soft function. The 
RG evolution shows that the color-correlated dPDF is 
Sudakov suppressed, in agreement with Ref. [37]. In ad- 
dition, we will also study the interference contributions 
such as shown in Fig. [5j which were first considered in 
Refs. [25j[26]. We find that all interference dPDFs are 
Sudakov suppressed. These conclusions greatly reduce 
the number of possible parton distributions that con- 
tribute to DPS at high energies. There is significant 



overlap of Sees. Ill and IV with the topics covered in 
a recent paper [26] which appeared while this work was 
in progress. Ref. [26] focuses on double parton scattering 
where the transverse momenta of the final state particles 
are measured. The evolution of the standard dPDFs, 
which we calculate in Sec. [V] differs from the evolution 
of transverse-momentum-dependent dPDFs obtained in 
Ref. [26] . This is similar to the situation in regular Drell- 
Yan, which is described by PDFs if the lepton transverse 
momentum is not measured (or integrated over), and by 
transverse- momentum-dependent PDFs if it is. The two 
are not simply related to each other, because they are 
renormalized differently. 

In Sec. [IT] we discuss the phenomenological aspects 
and implications of our results. The remainder of the 
paper contains the more technical aspects of our work. 
We present a systematic derivation of the factorization 
formula for the DPS contribution to double Drell-Yan 
production in Sec.|III[ after reviewing the steps that lead 



II. PHENOMENOLOGY 

In this section we give an overview of our work and its 
phenomenological implications. The technical details can 
be found in the remainder of the paper. As mentioned 
in the introduction, spin and color correlations can ap- 
pear for DPS. Already in SPS there are different spin 
dependent quark PDFs. For example, the polarized dis- 
tribution function Aq(x) measures the number of right 
minus the left-handed quarks with momentum fraction 
x, i.e. the longitudinal polarization. It contributes to 
the polarized structure function g±(x) that is measured 
in polarized deep-inelastic scattering. The transversity 
distribution hi(x) measures the transverse quark polar- 
ization Sq(x). It is a chiral-odd distribution, and does 
not contribute to polarized deep-inelastic scattering, but 
does contribute to polarized Drell-Yan with transversely 
polarized beams [38[ |39] . 

The quark spin in a proton is correlated with the 
hadron spin, so that both Aq(x) and Sq(x) vanish for un- 
polarized proton targets. By contrast, in DPS the spins 
of the two partons can be correlated with each other, and 
so nontrivial spin structures exist even if the proton is 
unpolarized. In addition to the usual unpolarized distri- 
bution F qq , there are two other distributions F^qAq, and 
FdqSq, which measure longitudinal and transverse spin 
correlations between two partons. F§ q § q is chiral odd in 
each 5q, but has overall chirality zero. 

In addition, dPDFs can also have color correlations 
which have no analog in SPS. The regular quark PDF 
is (schematically) the hadron matrix of the quark bilin- 
ear qq. For dPDFs, there are two possible color struc- 
tures qq qq and qT A q qT A q. Both objects are overall color 
singlets, but they give information on parton color cor- 
relations in the hadron target. We will refer to 1 <S> 1 
and T A <g> T A as the color-summed and color-correlated 
dPDFs. The two operators can also be written as lin- 
ear combinations of qqqq and q a q(3qPq a , which can be 
thought of as color-direct and color-exchange contribu- 
tions. 

Two quarks can be in a 6 or 3 color representation, 

which contribute to two different dPDFs Fqq and F qq \ 
measuring the color 6 and 3 diquark distributions. In 
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terms of the qqqq and qT A qqT A q basis, 
Ft q = 6F$+3F$, 

F T = 4C_p(6) _ 

qq \ qq qq J i 



(7) 



1) and F qq is 



where F qq is the color-summed dPDF (1 
the color-correlated dPDF (T A <8> T A ). (The factor of 4 
on the second line is due to an arbitrary normalization.) 
Since F qq and F qq are proportional to the probabilities to 
find diquarks in a color 6 or 3 and are positive, Eq. ( 7| im- 
plies -42^/3 < < 2F^/3. Similarly the qq dPDFs 
F^, F qq measure the singlet and octet qq distributions. 
The complete classification of dPDFs in term s of spin 

" FL 



and color structures is discussed in Sec. 



IV A 



and 



F qq evolve differently with energy, so the color correla- 
tions are energy dependent. 

In addition to spin and color correlations, there are 
also interference contributions to the cross section, an 
example of which is shown in Fig. [5] The interpretation 
of this as an interference contribution becomes clear in 
the context of a non-relativistic quark model, as discussed 



in Sec. [iVDl There we find that the regular dPDF F qt 



and the interference dPDF I qq - are (roughly) given by 



F 1 - 

qq 



\<t>qq(kl,k 2 )\ 2 



I 1 - 

qq 



^qqih.hy^qqikl.h 



(8) 

where <fi is the qq wave function and k\ and k 2 are the 
momenta of the quark and antiquark. The interference 
dPDFs do not have a nice probabilistic interpretation. 
Whereas regular dPDFs F(xi : x 2: z±) are positive (or 
real, in the case of spin and color correlations), the inter- 
ference dPDF does not even have to be real, as we argue 
in Sec. IV C| Only the contribution of the interference 
dPDFs to the cross section is real. In Sec. HVCl we also 
discuss the properties of dPDFs under discrete symme- 
tries. For example, charge conjugation invariance leads 
to F qq /p = F—/p, where the spin and color structures 
match up on both sides of this equation. 

To see how spin and color correlations and interference 
effects contribute to the cro ss section, requires a refine- 
ment of Eq. ([TJ. In Sec. Ill we systematically derive the 
formula for the DPS cross section for double Drell-Yan 
production, which results in the following leading-order 
factorization theorem 



da DPS 



( 



47ra 2 Qi 



dqfdY 1 dq 2 dY 2 V 3N c s 

2N r 



) ~J2J2 [ d2z ^l [i F qq F m + F AqAq F AqAq) + ( F qq F ^q + F AqAq F AqAq)] 
y H\H 2 J I 



■ * iy c r/ pT j?T ■ j?T W T \ , (W T W T , W T j?T \1 qTT 
1" \\ r qq r qq ~r~ r AqAq r AqAqJ "~ \ r qq r qq "T r AqAq r AqAq) \° 



2 

Nc 
' 2 

Nc_ 
~C F 



(4 



I AqAq) (Iqq 



AqAq) + IsqSqlsqSq 



sV 



(l£q + I AqAq) (Iqq + ^AqAq) + iJqSq^SqSq + 0- T) S] 



? T1 



(A 



L AqAq) (Iqq 



I AqAq) + ^ 



T rT 
dqdq 1 5q5q 



(9) 



Here qf and Y{ are the total invariant mass and rapid- 
ity of each lepton pair, and Q q is the quark charge. We 
suppressed the arguments of all functions in Eq. Q for 
brevity. The arguments of the dPDFs are the momen- 
tum fractions and transverse separation z_l, just as in 
Eq. 0. In addition to the dPDFs, the last four lines 
involve soft functions 5, describing the effects of soft ra- 
diation. The soft functions only depend on the large z± 
separation, since soft radiation does not resolve the short 
distances associated with the momentum fractions. The 
second line of Eq. ([9| contains the contribution from color 
correlations and the third through fifth line contains the 
interference contributions. Eq. ([9| agrees with the terms 
in the factorization theorem presented in Ref. [26] that 
were explicitly written out. 

The position space dPDF F(z±) is dimension two, and 
of order Aq CD . The momentum space dPDF F(r±) is 
dimensionless, and of order unity. The dPDF terms in 



Eq. 

- J d 2 z ± F(z ± )F(z ± ) (10) 

are of dimension two, and order Aq CD , since the z± in- 
tegral produces a 1/Aq CD . This shows that the dPDF 
cross-section is Aq CD /Q 2 suppressed. 

Clearly the large number of functions that appear in 
Eq. ([9| is worrisome for the prospect of measuring them 
in experiments. The good news is that color-correlations 
and interference contributions are Sudakov suppressed at 
high energies, leaving only the first line of this equation, 
which is one of the main conclusions of our paper. The 
Fs q sq dPDF does not enter in the leading order expression 
in Eq. 

Intuitively, the Sudakov suppression of color-correlated 
and interference dPDFs can be understood as a conse- 
quence of long range (zj_ ~ 1/Aqcd) color correlations. 
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FIG. 6. Color flow in the color-summed dPDF Fq q and color- 
correlated dPDF Fq q . The vertical separation of the vertices 
® is zi. 




FIG. 7. Color flow in the interference dPDFs Iqq and ijq. 
Both have long-range color correlations. The vertical separa- 
tion of the vertices is zi. 



This is illustrated in Fig.[6j where the color flow is shown 
for both the color-summed and color-correlated dPDF. In 
Fig. [7] we show the color flow for the interference dPDFs 
I 1 and 7 T , which both involve long range color correla- 
tions. While color is conserved in the hadron matrix ele- 
ment, the color-correlated and interference dPDFs move 
color a distance of order 1/Aqcd within the hadron. 

More formally, the Sudakov suppression follows from 
our study of the anomalous dimensions of the dPDFs 
and soft functions in Sec. El In the remainder of this 
section, we summarize our results for the RGE evolution 
and Sudakov suppression of F T . Similar results hold for 
the interference dPDFs, and are presented in Sec. [V| 

At this point we need to briefly discuss rapidity diver- 
gences, which we encounter in our calculations. These 
divergences cancel between the soft function and dPDFs 



and are thus absent in the cross-section. However, they 
need to be regulated and there is a corresponding series 
of large (single) logarithms that needs to be summed for 
reliable predictions. We achieve this using the recently 
introduced rapidity renormalization group [35 , 36 , whose 
workings are similar to that of dimensional regularization 
for UV divergences. Just as 1/e UV divergences lead to 
a fx anomalous dimension, 1/rj rapidity divergences lead 
to a v anomalous dimension. The new renormalization 
scale v is then used to sum the rapidity logarithms. (An 
introduction to the rapidity renormalization group and 
sample calculations involving the rapidity regulator are 



provided in Sees. VA and VC ) 

The anomalous dimensions for F T and S TT are calcu- 



lated in Sec. [VB] and Sec. |VC| and are given by 



F 



F 



xS(l — xi) 
a s (n)C A 



C F ~ \c^j P qq ( Xl ) + C A (\n^ + |) 



J(l-x 2 )^ 2 )(r ± ) + (1^2), 



1 



n 2 { rl/^) + 5{1 - Xl)5{1 - X2) 



7? 



20MCA U §{2) 

7T V 1 

tt 2a s (fi)C A 1 1 

2 (Jl /,,2\ ' 



tt 2 /i 2 (r 2 L //i 2 ) + 



(11) 



written in terms of the usual (one-dimensional) plus dis- 
tributions 



— = lim 



0(u - o 



u 



S(u-0 ln£ 



(12) 



In Eq. (11) we explicit showed the dependence on all 
variables, where a c)(l — x\) or means that the 

evolution does not affect the x\ or dependence, re- 
spectively. From Eq. ( pTj ) we read off that the natural 
scales for F T are (/iF,z/p) ~ (|rj_ | , p — ) ~ (A, Q) and for 
S TT are (jas^s) °° (A, A). A is a scale of order Aqcd 
which we take to be 1.4 GeV in our numerical analysis. 
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FIG. 9. Comparing the x-dependent part of the evolution 
of F 1 and F T assuming a common initial PDF. The three 
curves correspond to the initial PDF at \i — 1.4 GeV, and the 
evolution of F 1 and F T to fi = 100 GeV. 




Q 



FIG. 10. The Sudakov suppression factor with and without 
rapidity resummation, running from A = 1.4 GeV to Q. The 
evolution kernels U u and U u are discussed in the text. 



By evaluating F T and S TT at these scales, and running 
them to a common scale using the \i and v RGE, the 
large logarithms in the cross-section are summed. The 
natural scales and our running strategy are summarized 
in Fig. [8] As we mentioned, the rapidity divergences, and 
thus the corresponding v evolution, must cancel between 
the dPDFs and the soft function. Indeed, 



if + ^f T *(l - xi)5(l - x 2 ) 



0. 



(13) 



which provides a consistency check on our calculations. 

We will now show that the contribution of the color- 
correlated dPDF to the cross section is Sudakov sup- 
pressed. For both F T and 5 TT , \i ~ A is their natural 
scale, from we simultaneously evolve them to the hard 
scale /jl = Q. Combining the /i-evolution of the dPDFs 
and the soft function, we find 



7, 



1 gTT 



M ■ 2 

7T 



-7~ 5(1 -xi)5(l -x 2 ) 



CF- l -C A )p qq (x 1 ) + C A (\ii^ 



)S(1-X 1 ) 
(14) 



for each dPDF. The \i anomalous dimension has an in- 
teresting structure. The ^-dependent piece is the usual 
splitting function, but with a modified color factor of 
— 1/6 rather than 4/3, so the x-dependent part of the 
anomalous dimension can be interpreted as a slow "re- 
verse" evolution due to the change in magnitude and sign 
of the color factor. The effect of this reverse evolution is 
evaluated using; HOPPET g0] and shown in Fig. [9] for 
a sample initial PDF. It would be interesting to see this 
effect experimentally. 

There is also the second term in Eq. ( [14] ) , which is x- 
independent and thus does not change the shape of the 
dPDF. Its evolution (combining both dPDFs) from a low 
scale A to the hard scale Q ~ Qi ~ Q2 is given by an 



overall multiplicative factor 



17 M (A,Q) =exp 



In —J , 



2tt 



A 2 , 



(15) 



at leading-logarithmic accuracy. Color correlations are 
thus Sudakov suppressed and can be neglected for Q ^> 
A. Eq. (15) agrees with Ref. [37], which arrived at this 
conclusion by studying the color factors in the real and 
virtual Sudakov form factor. Ref. [37] did not attempt 
to factorize the cross section and does not discuss the 
interference case, which we find is also Sudakov sup- 
pressed. The effect of the Sudakov suppression is shown 
in Fig. [lO] at next-to-leading-logarithmic accuracy and 
taking the running of a s into account [which was ne- 
glected in Eq. (15)]. 

We also need to perform the ^-evolution to sum the 
rapidity logarithms. In Sec. |VD| we calculate the evolu- 
tion kernel U v . We estimate the effect of the ^-evolution 
on the cross section by (schematically) 

a~ J d 2 p ± d 2 q ± d 2 s ± F T (p ± ,^ F )F T (q ± ,^ F ) 
xU u (s ± ,v F ,v s )S TT (-p 1 _ - q_L - s_l,z/s) 

r\p±\<A 

- J d 2 P± ^(p±,Q,A), (16) 

Here we assumed that F T and S TT have a width of order 
A at their natural scales vp ~ Q, v$ ~ A. In Fig. 10 we 
take A = 1.4 GeV, and normalize our above estimate for 
the effect of the v evolution such that it is 1 when Q = A. 
The resummation of the rapidity logarithms increases the 
cross section, as shown in Fig. [lO] However, it is only 
a single-logarithmic series and thus has a smaller effect 
than the Sudakov suppression. At high energies the first 
line in Eq. ([9| thus dominates. 
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FIG. 11. Drell-Yan production, pi + p 2 ->> £ + £ 



III. FACTORIZATION THEOREM AT 
LEADING ORDER 

This section discusses the factorization theorem for 
DPS. We start by reviewing the usual derivation of the 
factorization theorem for single Drell-Yan, and then re- 
peat the analysis for double Drell-Yan. Other derivations 
of factorization for double parton scattering have been 
presented in Ref. [26] and for scalar partons in Ref. [5]. 



A. Single Drell-Yan 



The single Drell-Yan process is shown in Fig. 11 The 
Drell-Yan cross-section is 



da 



1 



2E P1 2E P2 v re \ 



d 3 h 



d 3 k 2 



(27r)32£ fcl (27r)32£ fc 



(27T) 4 d^(q-k 1 -k 2 )S^(q + 



PX ~ Pi ~ P2 



(17) 



where the sum is over hadronic final states X. The mo- 



menta are shown in Fig. 11 The amplitude for leptons 
produced through a virtual photon is 



A 



ieufaWvfa) iQ q (X\r\pp) , (18) 



where Q q is the electric charge of the quark. Performing 
the sum over the lepton spins 



spins 



4(k{k v 2 +k»k v x - g^k^) , 



(19) 



and the leptonic phase-space integral 
d 3 /ci d 3 fc ; 



{2n) 3 2E kl (2ir) 3 2E k 
x4(k^k^ + k^k^ - g^k 1 -k 2 ) 



(27T) 4 5^(q-k 1 -k 2 ) 



= ^ iff - fan 



(20) 



leads to 



(21) 



x (PP\4\X) (X\J v \pp) 1 (^f - g^) . 
The hadronic amplitude is shown in Fig. 12 and is 
(2^ 4 \q + p x - Pl -p 2 )M x 



(X\J»(z)\p lP2 ) 



= I d'z e iq - z e- iz < Pl+P2 - px) (A"|J"(0)lpip2> 
= (2n) 4 5^(q +p x - Pl - p 2 ) {X\J»(0)\ Pl p 2 ) . (22) 
The hadronic tensor can thus be rewritten as 

iP" = J2(^) 4S(4) (Q+Px -Pi -p 2 )M^M\ 



x 



= d 4 ze- iz - q (pp\J^(z)J v (0)\pp) , 



(23) 



such that 

4na 2 Q 



da 



3s 



f d 4 q 1 ( q»q v \ 



(24) 



The hadronic matrix element separates into two PDFs. 
It will be convenient to use light-cone coordinates where 
Pi is in the n direction and p 2 is in the n direction, such 
that Pi = n • pi and p~^ = n • p 2 are large. Explicitly, 
n M = (1,0,0,1) and n M = (1,0,0,-1) for beams along 
the third spatial direction. We find it convenient to use 
the same indices for both spin and color, where is 
Sab for the color indices, to reduce the number of indices. 
Since we do not observe the transverse momentum of the 
lepton pair, we can integrate over qx, 



(27T) 2 



cd 



= lablcd j 



dV 

(2tt)2 

dz + dz~d 2 zx 



d 4 z e -- + ? -/2- iz - g +/2 +iqi .z, fapi^WMzWMMVM 

i,+ 9 -/2-i,- 9 + / 2(5 (2) (z±) ^21^(^)^(^)^(0)^(0)1^2) . (25) 



At this point, one can contract the fields with the states in different ways. The total momentum at points and z 
have large — and + components of momentum. Thus one of the fields at z and must be contracted with pi and the 
other with p 2 . Furthermore, quark number is conserved, so one cannot contract two -0 fields with pi, etc. There are 
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two possible contractions which remain. The first one is 



= i^fcd I ^^e-- + ^/ 2 --"« + / 2 (p 1 |^Jz + ,0,0 ± )^(0)|p 1 )(p 2 |^(0,z-,0 ± )|p 2 )^ c (0) 



% lb % d >***f,(^)h(^) = + ^ - 2 ^<^Hk"> ■ (26) 



To derive the third line we note that z ± ~ 1/<2 T ~ 1/Q? 
whereas the dependence on z~ (z + ) of the matrix element 
of pi (P2) is slowly varying and may be set to zero. In 
the fourth line we have used the definition of the PDFs, 



47T 



fba 



^"/ 2 ( P |[^(^)W(^)] a [W t (0)^(0)] 6 |p) 

»£)'•£)-«(-£)*(-£): 



(27) 

The Wilson line goes to infinity and ensures gauge 

invariance, but was not written explicitly in Eq. (26) to 
avoid cumbersome notation. In momentum space the 
Wilson line is given by 



w n = ex p 



-9 



(28) 



where n denotes the direction of the energetic radiation. 



Combining Eq. (27) with Eq. (24) yields the familiar 
result 

da - da (29) 



dq 2 dY dq° dq s 



4ira 2 Q 2 



■ fq(Xl)fq(x 2 )] 



3N c q 2 s 

where we included the additional contribution to Eq. ( [25] ) 
from the second contraction with ijj ijj, which ex- 
changes the quark and antiquark PDFs. Here Xi are the 
momentum fractions and Y is the total rapidity of the 
lepton pair, 



Xi 



Pi 



n 2 

q e Y . 
s 



x 2 



pi 



In 2 

lq -e- Y . 
s 



(30) 



B. Double Drell-Yan 



The leptonic amplitudes for double Drell-Yan are iden- 
tical to that of single Drell-Yan, so the cross section is 
given by 



da = 2s 



3s J J (2tt) 4 (2tt 
(QiQi iuA 1 (Q2Q 

2 I „2 y ) ri 2 \ ri 2 



(2^) 4 (2^) 4 H ^ 



qf V qf J q$ \ q$ 



where the factor 2s ensures that we only count the flux 
factor once. For simplicity we will for now assume iden- 
tical quark flavors. The hadronic amplitude is shown in 



Fig. 13 and equals 



(27r) 4 5( 4 ) (q, + q2 +p x - Pl -p 2 ) Mf (32) 
= Jd i z 1 d 4 z 2 j*i-*i+*>-*2 (X\T{J?(z 1 )J 2 3 (z 2 )}\p 1 p 2 ), 



and the corresponding hadronic tensor is 





FIG. 12. Hadronic amplitude for single Drell-Yan. 



FIG. 13. Hadronic amplitude for double Drell-Yan. 
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H ^ap = ^2(2n) i S^(q 1 + q2+px - pi - p 2 )M x a *M x 

d 4 Z! d 4 z 2 d 4 z 3 e -i9i-*i-i»-* a +i*i-*s ( Pl p 2 \T{J^(z 1 )J^(z 2 )}T{J^z 3 )J^0)}\ Pl p 2 ) . (33) 



We now follow similar steps as for single Drell-Yan in order to factor the hadronic matrix element into dPDFs. The 
transverse momenta of the lepton pairs are not measured, so we can integrate over qi±,q2_L, 

/ f^LfS^L H ^ = I / d + d - d + d - d ^+ d - d 2 

J (2tt) 2 (2tt) 2 8 J 1 1 2 2 3 3 

(p 1 p 2 |T{J( lt (z 1) z ± )J 2 at (z 2) O ± )}T{Jr(z 3 ,z ± )j|(0)}|p 1 p 2 ) . (34) 

Note that one transverse position integral remains, in contrast to single Drell-Yan. This integral would remain even 
if one measures the transverse momenta of the lepton pairs. This can be seen from Fig. [2] the momentum r_\_ does 
not affect the transverse momentum of the final states, is arbitrary, and so is integrated over. The remaining d 2 z± 
integral in Eq. (34) is the position space version of this integral. It is important to note that is not an observable 
quantity [51] . 

The d 2 z± integral means that the correlation function of four currents, Eq. (34), is not local in z±. This allows 
soft radiation to contribute, since soft radiation can resolve the large distance z±. Soft radiation cannot resolve the 
short-distance scales zf ~ 1/Q. The soft radiation emitted by a quark field exponentiates into an eikonal (soft) 
Wilson line that can be factored out [4TJ [42] 

i/; f a (x) = S n , a > a {x)ip a {x) , (35) 

where the soft Wilson line is in the direction n of the quark and a, a' are color indices. 

There are various ways of joining the fields to the incoming particles in the matrix elements, which allow for large 
+ and — components of momentum flowing into the four current vertices. For example, 

f d 2 q i± f d^q2± ^ 

(36) 

= lla>b>~fc>d>le>f>lg>h> yW^rd4d2 2 -d4d2 3 -dVe- i9 i +Z i~ /2 - i9 i~ Z i +/2 e _i ^ ' Z * /2 - iq * ' Z l 12 e ^3-/2+i 9 r4/2 

X <p! \T {i> a ( Zl , z ± )^ c (z 2 , Ox) } T {V/ (z 3 , z ± )i> h (0) } \ Pl ) 
x <p2\T{tl> b (zu*±)M*2, 0±)} T {^ e (z 3 , Zx)^(0)} \ P2 ) 

x J dzfdz+dz+e-^^/ 2 e-^ z ^ 2 e^^/ 2 (pi\T {?/; a (zf , 0, z±)ip c (z 2 , 0, Ox)} T {ipf(z 3 ,0, z±)iph{ty} \pi) 
x j dz^dz 2 dz^e-^^l 2 e-^ z -l 2 e^^> 2 (p 2 \T {^(0, zf, z ± )^(0,zj, Ox)} T {^ e (0, z£ , zx)V^(O)} \p 2 ) 
x (0|^,aa'(0, 0, z ± )SnM^, 0, z ± )4 cc , (0)5^(0)4^, (0, 0, z ± )S nJ , f (0, 0, z ± )Sl gg ,(0)S n , h , h (0)\0) . 

Here we once again use that zf ~ \/Q varies rapidly, whereas in the matrix element of p\ (^2) the dependence on 
z~ (zf) is slow and can be set to zero. The dependence of the soft radiation is slow in both zf and z~ and so only 
the dependence on z± remains. The last line in Eq. (36) is the soft function S (we omit a subscript qqqq since we 
do not consider soft functions with gluons in this paper)7 The soft function has a lot of indices but also has a lot of 
symmetry. When we contract indices below, only four independent soft functions will appear. We identify the dPDFs 
F qq and F m in Eq. Q by 

F qq { q ^q 2 ^p- 1 ) ac}h = -i^p- 1 J ^^L^L e -^4/* e -^4/2 e ^4/i 

( Pl \T{^ a (z+, 0,zx)? c (4>0,Ox)}T{^ / (4,0,zx)V/ ( (0)} \ Pl ) , 

F m ( q t,q},z ± , P l) bdeg = -Wn P l J ^^^e-^r/2 e -i^ 2 /2 e i^ 3 -/2 

(p 2 |r{^(0^ 1 -,zx)^(0,z 2 -,Ox)}r{^ e (0,z 3 -,zx)^(0)}|p 2 ) . (37) 
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The factors p± and p 2 are included to make the dPDFs boost invariant, indicating that each dPDF comes with a 
Aqcd/Q factor in the cross-section. These dPDFs are dimension two objects and of order Aq CD (including the p^ 
factor). The overall minus signs in Eq. (|37|) are due to the ordering of fermion fields. Translation invariance gives 



(38) 



Fqqiq.1 1Q2 , z ±,Pl)acfh = Fqq{q.2 \ Ql "j ~ z ^,Pl)cahf 1 
Fqqiq.1 ,0.2 , z ±,Pl)bdeg = Fqqifh 1 Ql > _Z -L j Pi )dbge • 

The other possible ways of connecting the fields with the states in Eq. ( [34] ) that contribute are 

1 J 4tt 4tt 47r 
(pi |T {^ a (*+ 0, z ± )^ d (4, 0, ± )} T {^/(4, 0, z ± )^(0)} |pi> 



^Wfe ^2 5 Z ±,Pl W<? 
F qq{Ql,q2i Z ±iPl)bceh = -16tT^ 



y dz+ d4 

./ 47T 47T 



d4 

47T 



2+ /2 -i<? 2 Z+/2 ig 1 z+/2 



hq^l ^2 ^ Z ^P\ )adeh = lfa P 1 



tqqiQl ^ Z ^P\ )bcfg = ^ Pi 



(pi \T ty b (zt, 0, zx)V c (4, 0. ± )} T {^ c (4> °> Z ±)^(0)} |P1> 

dz^dz^dz^ e - iq - z + /2 e -iq~z+/2 jq^ z+ /2 
47T 47T 47T 

(P! |T {^ a (4, 0, z ± )r/; d (z+, 0, Ox)} T {^ c (4, 0, z ± )^(0)} |pi> 

dzj" dz£ iq7z + /2 iq - z +/2 iq-zt/2 



4tt 4tt 4tt 



<pi |T {^(z+, 0, zx)^ c (4» 0, 0j_)} T 0, z ± )^(0)} bi) 



(39) 



is the interference double PDF, corresponding to the contribution shown in Fig. [5] The origin of this name will 
become clearer in Sec. IIVDI Translation invariance leads to the relations 

F qq(<ll,<l2i Z -LiPl)adfg = Fqq{0.2 i Ql > ~ Z ±- iPl)dagf , 

I qq(QliQ2i Z -LiPl)adeh = Iqq (fe ,Ql , ~ z ± , Pl)dahe • (40) 

The soft functions that correspond to these other contractions can be expressed in terms of the soft function S. For 
FqqFqq and Iqqlqq respectively, we have 

(°l^,aa'( Z ^)£n,6'6^ = S aa 'bb'gg' hh'ee'ff cc> dd> , 

<0|< aa ,(zx)^ 6 (zx^^^^ 'bb' ee' f f cc' dd' gg' hh' • 

(41) 

Combining all these ingredients, we find that the cross section is given by 



do- 



DPS 



dq\ dYi dg| dY 2 



X {.Fqq^cfhFqq^deg^aa'bb'cc'dd'ee'ff'gg'hh' H~ F qq,adf gFqq,b 



bceh^aa'bb' gg' hh' ee' f f cc' dd' 



(1) 



r(2) 



qq^adeh 1 qq,bcf g^aa'bb' ee' f f' cc' dd' gg' hh' 



+ (9 ?)] , 



(42) 



where the first dPDF is for the first hadron, and the second for the se cond h adron. This expression can be simplified 



IV A 



(This is analogous to e.g. (pi |V y0b|pi 



using the color and spin decompositions of the dPDFs, discussed in Sec. 
dab by color invariance, which we encountered in single parton scattering.) Using the decomposition in Sec. 
find the result presented earlier, 



IV A 



we 



da DPS 



dq( dYi dqi dY 2 



= ( 



4ira 2 Q 2 „\2 1 



V 3N c s 
2N C 



) ~J2J2 [ d 2 Z± \ [(FqqFfq + F AqAq F AqAq) + ( F qq F ^q + F AqAq F AqAq)] S ^ 
' H\H2 J \ 



c F 



[(FqqFq^q + FAqAqFAqAq) + ^qq^lq + F AqAq F AqAq )]S 



1 

2 

Nc 
' 2 
Nc 

~C F l 



(I^q + ^AqAq)(^lq + ^ 



L AqAqJ 



rl_ rl 

5qSq Sqdq 



s} 1 



(l£ q + lAqAq)(Iq~q + ^AqAq) + ^Iq5q^5q5q + (1 ^ ^) 



(/^ + lAqAq)(Iqq + ^AqAq) + ^IqSq^IqSq + (Q ^ Q) 



(43) 
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The soft functions are denned below in Eqs. (44)- (46) and are normalized to be 1 at tree- level (except for Sj 1 which 
vanishes at tree-level). We suppressed the arguments of the functions for brevity. The arguments of the dPDFs are the 
transverse separation z± and the momentum fractions, given by equations analogous to Eq. J30| ). The soft functions 
only depend on z±. Eq. (43) reduces to the familiar result in the literature given in Eq. dlj), when spin and color 
correlations as well as the interference dPDF are ignored. A rigorous proof of factorization requires one to show that 
Glauber gluons do not contribute, which we assume here. 

The factorized form Eq. (43) contains soft functions from the soft Wilson lines in Eq. (35). The soft function for 
the color-summed dPDF is 



S 11 — "T72 dfadhc $be$dg $a'b' $c' d' $e' f $ g' h' S qqqq,aa'bb' cc' dd' ee' f f' gg' hh' 
^ c 

= ^ <0|tr[Sj(z ± )SW(z±)4(^)5n(z±)]tr[St(0)SW(0)4(0)5„(0)]|0> = 1. 



(44) 



using unitarity of the soft Wilson line. The color-summed soft function S 11 is trivial and receives no QCD corrections. 
The soft function for the color-correlated dPDF is 



C F N C /! 
2 



C F N C 



T fa T hc T be T dg $a'b'$c' d'$e' f'$g'h' Sqqqq,aa'bb> cc> dd> ee' f f> gg> hh> 

<0|tr[T^t(z ± )S^ 



(45) 



The color-correlated soft function is nontrivial, and has been normalized to unity at tree-level, S TT — 1 + (D(a s ). The 
soft functions for the interference terms are given by 

S V = ^2 <0|tr[4(z ± )Si r (z ± )4(0)5 n (0)]tr[St(z ± )SW(z ± )4(0)5 n (0)]|0> , 



si 



<T1 



sj 



TT 



C F N* 



C F N C 



i|tr[T A 4(z ± )SW(z±)4(0)5„(0)]tr[T A 4(z ± )S 5r (z ± )4(0)5 n (0)]|0> , 
(0|tr[T A 1 5t(z ± )^(z ± )T B 4(0),5„(0)]tr[T' 4 4(z ± ) > S 7r (z ± )T B 4(0)5„(0)]|0) 



(46) 



Defining 

W = 4(z x )5 n (z x )4(0)SW(0) 

y = 4(z x )sw(z x )4(o)5 n (o) . 

and using the identity 



1 



ArpA 



TIT- 



b x cd 



^adScb — -^J^dabdcd , 



they can be written as 
O'J 

1 



— 1— (0\ti[T A WT A W^\0) 
C F N C 



s} 1 
sp 



2C F N t 
^2 <0|tr[^] 2 |0> , 



(OltrlWjtrlW^IO) 



2C F N C 



C F N* 
1 



MT A y] ti[T A y}\0) , 



<0|tr[}f|0> 



1 



C F NI 

1 (o\tv[T A yT A y}\o) 



C F N? 



{o\tv[y 2 



C F N C 

1 



C F N* 



<\tv[T A y} tx[T A y]\o) 



(47) 
(48) 

(49) 



2C F Ni 



<0|trM 2 |0> 



|tr[y 2 ]|0) . 



None of these soft functions are trivial. S TT ', S} 1 and 
Sj T are normalized to 1 at tree-level. The soft function 
Sj 1 connects interference dPDFs with color structures 
T A (g) T A and 1(8)1, and only starts to contributes at 
order a s . The three interference soft functions are not 
independent, but satisfy 

Sf^Sf-Sr, (50) 

so only two independent interference soft functions exist. 

IV. DOUBLE PDF 



In this section, we classify the allowed color and spin 
structures for the dPDF. The F qq dPDF matrix element 
has the schematic form 

(P\^a^c^f^h\p) , (51) 

where the subscripts represent color and spin indices, 
with similar expressions for the other dPDFs given in 
Sec. IIIIBI The allowed color and Lorentz structures for 
these four-quark matrix elements are derived in this 
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section. The possibility of nontrivial color and spin 
structures for unpolarized dPDFs wa s first discussed in 
Ref. [27 j. Our decomposition in Sec. IV A is essentially 



ble PDFs of unpolarized hadrons, 



the same as in Ref. [26 , but the normalizations are 
slightly different. 



A. Classification 



The qq dPDF has two possible color structures, since 
the proton is a color singlet state, 



(52) 



These are used to decompose the double PDF into the 
color-summed and color-correlated dPDFs as 



(Fqq) ac f h 

(%) 

(Fqq) ad f g 

(Fqq) bceh 



Jj2 S faShcFqq + ~C^~^fa T hc F qq , 

bdeg ~~ Jp^eSdgFqq + (J p ^ J^be^dg Fqq > 

1 r r rpl | ^ A A T 
N 2°f a ° d 9 r qq^ q 1 fa- L dg r qq > 



F^c 

1 1 2 

m SbeShcF ™ + Cf~n, 



TfoThc Flfq J (53) 



and similarly for I g g and 7^. The double PDF with the 
T A (g) T A color structure can be interpreted as measuring 
diparton color correlations. This is clear from decompos- 
ing F qq into the 6 (symmetric tensor) and 3 (antisym- 
metric tensor) and F qq - into the 1 and 8 (adjoint) color 
structures: 



7?(6) 

qq 



7?(3) 



,(8) 



— — F 1 
Af 2 qq 

1 y c 

N 2 99 



1 


+ 1) 


N C (N C 


-1) 






N c qq 




2 




N C (N* 


-1) 



qq ' 



(54) 



( n a <7^ a 75 ) a / ( % CT^ 75 ) c/l 



(55) 

-(^7±75)a/(^7u75)ch, 



In SPS only the (analog of the) first spin structure con- 
tributes for unpolarized hadrons. The other two struc- 
tures enter only for longitudinally and transversely polar- 
ized protons, respectively. They contribute to the g\{x) 
polarized structure function in DIS, and the transversity 
distribution h\{x), respectively. In the dPDF we have 
two partons, and the last two structures can appear for 
unpolarized protons due to diparton spin correlations. 
The decomposition of the dPDF into the various spin 
structures is given by 



1 



( F qq)acfh = Jq [$fa$hc F Q<I + (fa) fa (fa)h C F A qA q] 



+ ^ a l^lb)fa(°^ + lb)hcF5q8q • 



(56) 



The numerical factors have been chosen to match with 
the conventional normalization of the PDFs. 

To clarify the spin structures, it is helpful to write 
out these double PDFs in terms of quark creation and 
annihilation operators 



F ai 



qq - (p\( a \R a iR + a\ L a 1L )(al R a 2 R + a\ L a 2L )\p) , 

FAqAq ~ <j)\(a[ R a 1R - a\ L a 1L )(a\ R a 2R - a\ L a 2L )\p) , 

FsqSq ~ (p|2(a{ i? ai L 4 I/ a 2 i? + a| I/ ai jR 4 i? a 2 L)b) , (57) 

where we have assumed that the proton is in the n M = 
(1,0,0,1) direction. F qq measures the joint probability 
to find two partons in the proton, F Aq A q measures the 
correlation between the longitudinal polarization of two 
partons in the proton, and F$ q 5 q measures the correlation 
between the transverse polarization of two partons in the 
proton. FsqSq does not enter the leading-order factoriza- 
tion formula Eq. ([9|, but can contribute at higher order 
in a*. 



B. Definitions 

Starting from the dPDFs with uncontracted indices 
in Eqs. (36) and (39), we now apply the spin and color 



decompositions of Sec. |IV A] This leads to a large number 
of dPDFs, most of which contribute to the cross section in 



Eq. ( 43 ) . We show a small sample of these for illustrative 



There are three spin structures that appear in the dou- purposes below: 



Fi(% > % >a± )=-167rpr 
K Pi Pi J 



dz^dz^dz^ e - iq - z + /2 e -iq~z + /2 jq^ z£ /2 
47T 47T 47T 



f 



( Pl \T\ V(z+0,z ± )|j o V(4.0,0±)^ JT{ifj a (4Az±)M0)}\Pi) , 



2h 
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2^(%,%,zl)=-167t Pi / — ^--e 



dz^dz^dz^_ -- z + /2 iq - z + /2 iq-z+/2 



Att Att Att 

<pi |T | [^(*+ 0, z ± ) |t a ] ^ 6 (4, 0, ± ) J T |^ a (4, 0, zx) [^(0)|t a J & j> \ Pl ) , 



_) = — 16^^! y 4 



^-p-^i 4 i 2 P ~ i( ii 4 i 2 P i( ii z i i 2 



Att Ait Ait 



dz£ dz£ dz4 
~A 

,+ 
a 

\Pl Pi / ./ 47T 



<pi|t{[#z+ o,z ± )^^^ , ^h) , 



47T 47T 



^g(%,%,zi) = 167r Pl 



(pi|t|^ 6 (4,0,z_l 

dz^ d^2~ dz ; + 



JT |Va(4.0, 



'3- P - i( ii z t/ 2 p-^ z i i 2 p i( i\ z i i 2 



Att Att Att 



(Pi\T 



0,z ± ) ^(4,0,0 ± )^- 7 5 



o |t|vo(4.o, 



z±) V>(o)-p75j 6 fbi) • 



(58) 



We have written the arguments as q± /p± and ->P2 •> 
since by boost invariance the dPDFs can only depend 
on these combinations. In addition to the two color and 
three spin structures, we also have different quark flavors, 
which have been suppressed. For example, one can have 
dPDFs F uu , F ud , F AuAd , etc. 



C. Discrete symmetries 

In this section we study the properties of the double 
PDFs under discrete symmetries, for which it will be con- 
venient to use the dPDFs with uncontracted indices. We 
start by considering the complex conjugate of the double 
PDF: 



F (a~ a' z, T)~Y *u--16nn7 f e ^4l 2 e ^4l 2 e ~^4/ 2 

fqqKQl ,<12 ^^V\ Jacfh ~ ^ Pl J ^ ^ ^ e e e 

x ( Pl \T { fiftOh ] h Mz£ , 0, zx ) 7 °] / } T { [7V(4,0,0 ± )] c [ 7 V(4,0,z±)]a}bi> 

= -16ttv7 [ d ^ d ^ d ^t e -^i4/ 2 e -^2 4/ 2 e ^i4/ 2 
1 J Att Att Att 

x ( Pl \T { f0(*+ 0, zx)7 ]/[^(4 5 0, 0±)7°]4 T { [ 7 V(4> 0, z±)] a [7V(0)] c } \pi) 

= F qq(Ql,Q2i Z ±iPl)acfh • (59) 

In the first step we interchange z\ «->• 2:3 and use momentum conservation to change the field in which the z£ coordinate 
appears. In the second step we use 



7°(#) T 7° 



t, 7 (^75) T 7°=^75, 7°(^75) T 7 = ±^75, (T A ) T ® (T A ) T = T A <g> T A , (60) 



to show that the spin and color structures are unaffected. We thus conclude that F qq is real. This is true for the 
other double PDFs, except for I qq and I qq which satisfy 

kq(<h \Q2 '^^PlT = Iqq(QliQ2i Z -LiPl) • ( 61 ) 

Thus the individual interference double PDFs are not necessarily real, but their contribution to the cross section is. 



This is not surprising, given their interpretation in the context of a quark model in Sec. IV D 
Under parity 

F qq ( qi ,€,*±,Pl)U = -16^r / m^" 4/2 ^ 4/2 ^ 4/2 (62) 

x(pf|T{W0,z^-z ± )7 O ]a[^4,O i )7 ] c }T{[ 7 V(0,4,-^)]/[7V(0)]/ l }lpf) • 
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Using equations analogous to Eq. (60), we find that the spin and color structures are unchanged (with n <H> n), 
yielding F qq for a proton in the n direction with z± — » — z±. This is not particularly useful, since the double PDF for 
a proton in the n and n direction are already related by a rotation. 
Under charge conjugation 



F qq/p(Ql ,Q2 ,Z±,Pl ) 



c 

acfh 



-167rp 1 



^4 ^4 ig- 2 +/2 p -i q -4/2 ig-4/2 



47T 47T 47T 



x (p, \T { [ 7 Vfy(*i + , 0, z ± )]^7VV(4, 0. 0±)7°]H T { [^(z 3 +, 0, z ± ) 7 V]J[^0) 7 V]£} l^i> 

F qq/P^ ^2^^ Pi) fhac • (63) 



relating the double PDF for the proton and the antiproton. Using equations analogous to Eq. (60), we find that the 
spin and color structures are unchanged, e.g. F AqAq ^ p — F AqAq -/p- 

I 

The matrix element is 

5 r , s 5 r ,, s ,(27r) 3 5( 3 \k' - k)(2^) 3 ^ 3 )(k' - l')(27r) 3 <J( 3 )(k - 1) 



D. Interpretation 

We will now provide an interpretation of the double 
PDF using a nonrelativistic quark model. Alternatively, 
the interpretation of dPDFs becomes clearer when they 
are written in terms of light-cone wave functions of the 
colliding hadrons [43 . As a warm up we start with the 
single PDF in Eq. ([27]). Since it is boost invariant, we 
can work in a frame where p is at rest, 

/d 3 k 
^3<A 9 (k, S )|k, S )®|p-k).(64) 

Here |k, s) is a quark with momentum k and spin/color 
5, and |p — k) represents everything else that makes up 
the total momentum of the proton. We have switched 
to a nonrelativistic normalization for the states. The 
normalization of the wave function <j) q is given by 



4-7 (27T 



^|0,(k,s) 



1. 



(65) 



Writing the fields in the PDF in terms of creation and 
annihilation operators, 



'«(f) 



2E n 



47T 



-iz + q~ /2 



x (k',s' 



d 3 l 

(2tt) 3 



31/ 



d a l 

(2tt) 3 



z+l'~/2 



> p — k'| ^'^ a ^ |k, a Op — k) 



x«(l',r')|«(l,r). 



2£^ 



leading to 



(67) 



/tr e " iz+9 " /2 /(0'^ (k)|2 ^" 72 



x u(k, s') — u(k, s) 

2E k ^ v 7 2 v 7 



p- 



■) • (68) 



This gives the single PDF its probabilistic interpretation. 
It is the probability to find a parton in the proton for a 
given value of x = q~ jp~ . 

We will now show that a similar probabilistic interpre- 
tation is possible for the double PDF. Boost invariance 
allows us to go to the rest frame again, where we write 
the wave function for two quarks in a proton as 



— » 



/2E P ^ fd 3 k, fd 3 k 2 



si,s 2 

x |ki, 81, k 2 , 5 2 , p - ki - k 2 ) . 



(69) 



We are assuming a single quark flavor for simplicity, and 
accordingly included an overall factor of 1 /2 in the above 
equation. The normalization of the wave function (p qq is 

\T.j^Sw?^''^ s ^ = ^ (70) 



(66) Following the same steps as for the single PDF, we find 



1 / N A f d 3 ki d 3 k 2 / 1 1 \ * , / 1 1 

F qq (x 1 ,x 2 ,r±) = 4^ J 7^)3 (2^)3 ^ V 1 + 2 rj -' Si; 2 2 T± ' S<2 ) ^ qq \ 1 2 r± ' 51 ' 2 + 2 T± ' S<2 



xS(xi—k 1 /p )5(x 2 -k 2 /p ). 



(71) 



16 



We see that corresponds to a transfer of transverse momentum between the two double PDFs (see Fig. [2|. In 
this form the double PDF does not have a probabilistic interpretation. However, if we first switch to a wave function 
that depends on k~ and (in which we absorb the Jacobian) and then Fourier transform to position space for all 
transverse variables, the dPDF does have a probabilistic interpretation, 



dk 1 d 2 kix dk 2 d 2 k 2 ± 
3 (2tt) 3 



(k t , k i± + ^r_L, s i; k 2 , k 2± - ^±^2) 



xcfrqq^ ,ki ± - ^r_L,si; k 2 , k 2± + ^rj_, s 2 )^i - K /p )S(x 2 -k 2 /p ) 

/ \ 2 r» 2 



(72) 



Likewise, for the other double PDFs we find 

d 3 ki d 3 k 2 ± 



( 1 1 \* / 1 1 \ 

r(ki + -r_L,si;k 2 - -rj_ 5 s 2 J 0</g(ki - ~r ± ,si;k 2 + ^'^J 



F^(x u x 2 ,r ± ) =4^|pp ^ 

X(5(xi - /cf/p")(5(x 2 - k 2 /p~) , 

1-1/ x , v- /"d 3 ki d 3 k 2 ^ / 1 ,1 \* / 1 ,1 

J^(xi,x 2 ,r_L) =42^/ 0<ra(k 2 - -r_L,s 2 ;ki + ^r^si J 0<rc(ki ^r^,^;^ + 2 r± ' 52 

x5(xi - fcf /p~)S(x 2 - k 2 /p~) . 



(73) 



By swapping the momenta, we can write I qq - as (jj^ q (j)qq, 
which makes it clear why it is called the interference dou- 
ble PDF. The interpretation of the various spin and color 
correlations was given in Eqs. (54) and (57). 



V. RENORMALIZATION GROUP EVOLUTION 

In this section, we compute the renormalization group 
evolution of the dPDFs and soft functions using a re- 
cently introduced rapidity regulator [35j [36] . 



A. Rapidity RGE 

In the calculation of the RG evolution of the dPDFs 
and soft functions, we will encounter so-called rapidity 
divergences. These arise because the collinear degrees of 
freedom (the dPDFs) and soft degrees of freedom (soft 
functions) are only separated in rapidity and not in in- 
variant mass. Typically, rapidity divergences arise in in- 
tegrals such as 



f Q dk^ __ f p dk^ f 

Jk ± k~ J k± k~ J p 

r°° dk~ | r 

Jk 1 k~ Jo 



dk~ 

Q dk2 
k- 



(74) 



In the factorized cross section this integral gets split into 
a contribution from the soft and from the collinear re- 
gion, as shown on the first line. On the second line 
we systematically expand in the power counting, where 
in the soft region p — >> 00 and in the collinear region 



p —> 0. Both the soft and collinear contributions have 
a rapidity divergence. Rapidity divergences only appear 
when one factorizes the cross section, and must cancel 
between the collinear and soft contributions. These di- 
vergences need to be regulated and the corresponding 
series of large (single) logarithms of k±/Q in the cross 
section need to be resummed for a reliable prediction. We 
achieve this using the recently developed rapidity renor- 
malization group [35j [36] , which was introduced in the 
framework of soft-collinear effective theory [42j [44H46] . 

Note that factorization already splits momentum inte- 
grals by invariant mass, 



r Q d A k _ r p d A k r Q 

Ja if "A ^ + i P 

r 00 d A k r 

"A Wo 



d 4 fc 

3 d 4 fc 
fc 4 



(75) 



where A is an infrared scale of order Aqcd- One may 
think of the first and second term as (roughly) corre- 
sponding to the PDF and the partonic cross section. In 
dimensional regularization the ultraviolet divergence of 
the first term cancels the IR divergence of the second 
term. The ultraviolet divergence leads to an anomalous 
dimension (for the PDF this would yield the usual split- 
ting functions), which can then be used to sum the large 
logarithms of Q/A. The problem with rapidity diver- 
gences is that they are not regulated by dimensional reg- 
ularization. 

The rapidity renormalization group modifies the 



collinear Wilson lines [see Eq. (28)] and the soft Wilson 
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FIG. 14. Real and virtual contributions to the PDF evolution. 
The Wilson line is denoted by a double line. 



lines, 



^ = 1-^-^(0). 



gw v^l 2 
id+(2id 3 )i/ 2 



(76) 



Here rj is the rapidity regulator, v a new renormalization 
scale and w = 1 is a bookkeeping parameter. The proper 
extension of Eq. (76) beyond one-loop, where you can 
have multiple emissions, is discussed in Refs. [35j [36] . 
Rapidity divergences arise as 1/rj poles in calculations 
and will lead to a ^-anomalous dimension, in complete 
analogy to 1/e poles and the \i anomalous dimension. 
The rapidity logarithms will be summed using the v- 
RGE [13136]. 

We will illustrate the use of the rapidity regulator by 
calculating two diagrams that arise for the single PDF. 
The PDF has no rapidity divergences, so the rapidity di- 
vergences cancel when you add these diagrams. However, 
these same diagrams appear with different color factors 
for the color-correlated dPDF. In that case the rapidity 
divergences will no longer cancel. We start by consider- 
ing the left graph in Fig. 14 Summing over the external 



polarization and introducing a gluon mass M to regulate 
the IR divergences, 



p — k ' 



Itr 



\- 9™ 2 n p ^T A t ifr-fl pt a 1 
(fc-) 1 ^ 2 (p-k) 2 yi p \ 



_ g 2 w 2 C F (H 2 e IE y / v \v x rd d ~ 2 kj_ 1 

~ 2^ V 4tt / Kp^) (l-xy+i J (2w) d - 2 k 2 ± +xm 2 

a s w 2 c F (fycY xl ~ e 

2tt y > \m 2 J \p-J (1-xY+v 

= ^ { [ - -r( £ )e- (4V - ' In -] 6(1 - x) + 1 —V + 0(e%°) ) . 
2tt { L r] \m z J e p J e (1 — x)+ J 



(77) 



To obtain the expanded expression on the last line, we 
used the distribution identity 



T, "TTT" = ~-S(l - x) + j- 1 r . (78) 

(1 -xY+v r] v ' (1 -x)+ y J 



[The definition of the (standard) plus distribution was 
given in Eq. (12).] 



Absorbing the divergences on the last line of Eq. ( 77 ) 



into the renormalization factor Z\a-> the contribution 
from this diagram to the anomalous dimension is 



dZ 



1A 



a s C F 



x v 
5(1 — x) In — 



11B 



2ig 2 w 2 C F 5(l 
g 2 w 2 C F 



5(1 - x) 



f d d k 
JW) 



dZ 



1A 



dv 



a s C F n 2 



x). 



(79) 



Here we used 



da s 
dfi 



= -2ea s +O(e ), 



dw 



'■w + 0(rf>), (80) 



which illustrates the role of w. Graph 1A also has a 
mirror image and thus contributes with a combinatorial 
weight of 2. 

Note that the anomalous dimension contains an ex- 
plicit dependence on p~ , indicating the breaking of boost 
invariance. This is required, since any rapidity regulator 
must break boost invariance [see e.g. Eq. (74)]. Of course 
the cross section is still boost invariance. 



The second graph in Fig. 14 is 



tr 



t gw 2 n p ^T A \(f - #) . 



2tt 
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a s w 2 Cp 
2^ 

a s w 2 CF 



*(i-.)r<*~(£,V(£)' 



2tt 



S(l-x)T(e)e 



\M 2 J \p- 

(j*y 



du ■ 



,1+77 



2tt 



\M 2 J 77(1 — 77) \p 

5(1 -x) -r(e)e e ^(^4V + -ln- 
r\ Vmv e p 



a s w 2 CF x fl 



0(6V) 



(81) 



In the second step we performed the /c + integral by con- 
tours. This leads to < k~ < p~ , since the poles are 
otherwise on the same side of the real axis. We then 
subsequently perform the and k~ integrals, and ex- 
pand in r] and e. The contribution of this diagram to the 
anomalous dimension is given by 

7iB,, = ^(ln- + lW-*) 

7T \ V J 



p~ 

a s C F 11 2 



(82) 



This graph also has a mirror image and so has a combina- 
torial weight of 2. As we anticipated, Jia,v + 7ib,v = 0. 
However, this will no longer be true for color-correlated 
and interference double PDFs. In those cases, we get 
essentially the same graphs as discussed here. The first 
graph still has color factor CU, but the second graph has 
a different color factor, e.g. Cp — Ca/2 for the color- 
correlated dPDF. The v anomalous dimension no longer 
vanishes for the dPDF, but cancels against the v anoma- 
lous dimension of the soft function. 



B. Double PDF 

In the previous section we gave a brief introduction 
to rapidity regulator and illustrated its use for two ex- 
plicit examples. Here we simply tabulate our results for 
the various diagrams. The diagrams involving a single 
quark line are shown in table [T| with separate columns 
listing the color factors for the 1 (g) 1 and T A <g> T A color 
structures. Since no transverse momentum can be trans- 
ferred these graphs are all proportional to S^ 2 \r±). The 
diagrams connecting both quark lines are shown in ta- 
ble [Hj In principle these diagrams allow for mixing be- 
tween color structures, but this cancels in the sum over 
diagrams. Here we used the color factors Cd and C±, 



c d 



N 2 



N 



N 



1 



4/V 2 



(83) 



and the standard plus distribution defined in Eq. (12). 



Adding up these diagrams, we obtain the anomalous 



d f dy 1 dy 2 2 , fxi x 2 

fi— F(x 1 ,x 2 ,y 1 _) = / d ki7 M — ,-,ri-kj 

d/i J 2/1 y2 \yi 2/2 



dimensions of the dPDFs. The corresponding RGE is 

' dj/i dy 2 

2/1 2/2 " v ?7i 2/2 

xF(2/i,2/2,k_L 

d 1?( , fdy 1 dy 2 2 /xj x 2 

v—F(x 1 ,x 2 ,r±) = / d k_L7M 

dv J 2/1 2/2 V 

x%i,l/2,ki) 



^2/1 2/2 



,r_L-k 





(84) 

where we have chosen not to put any factors of 2tt in the 
convolution integral. It is convenient to write our one- 
loop results for the dPDF anomalous dimensions using 
% t1/ defined by 

7m (*i, * 2 , r ± ) = ^ (r ± ) [%( Xl )S(l - x 2 ) 

+5(1 - x 1 )%(x 2 )] , 
jv(xi,X2,r±) = %(r ± )S(l-x 1 )5(l-x 2 ) . (85) 

7T 

For the color-summed dPDF F 1 , the total anomalous 
dimensions are 



7f {*)=C. 



(1 

C F P qq (x) , 



^. +( i-x)c + |*(i-x); 



ll (r±)=0. 



(86) 



which is the usual evolution for each of the individual 
quarks. It is understood that the splitting function is 
modified, depending on the spin structure. The splitting 
function for the longitudinal polarization is the same as 
for the unpolarized case, PAqAq(x) = P qq (x). For the 
transverse polarized case, P5 q 5 q (x) — P qq (x) — (1 — x). 
For the color-correlated dPDF we then find 

1 



% (*) =[C F --C A P qq (x) 



+Cu(ln^ + -x), 



F (rx) 



3^ 

P~ 4> 

C A 1 1 

7T \± 2 (r2_//i 2 ) + ' 



(87) 



Note that the gluon mass M drops out in the sum over 
diagrams and does not appear in the anomalous dimen- 
sions, as it should. A simple cross-check on these results 
is provided by 
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TABLE I. Double PDF renormalization from diagrams involving a single quark line. The columns show the graph, the color 
factors for the 101 and T A 0T A color structures, the ji and v anomalous dimension and the weight factor for each contribution. 
Only Graph II depends on the spin structure through £, where f = 1 for q, q, Aq, Aq and £ = for 5q, Sq. 
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T A 0T A 


7m 
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0(10 1) 
T A 0T A 


!(C d + CU)T A 0T A 
Ci 1 1 





1 \ f 2 £ (2) (r±)ln 
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1 \, 2 } a , +5 (2) (r±)ln j& 
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TABLE II. Double PDF renormalization from diagrams involving both quark lines. The columns show the graph, the color 
factors for the 1 1 and T A T A color structures, the \i and v anomalous dimension and the weight factor. The graphs do 
not depend on spin. The diagrams are not diagonal in color, and the anomalous dimensions are weighted by the sum of both 
color factors for a given graph. 
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which can be verified using 



d 1 



= -25{r\) = -27rS^(r ± ). (89) 



From the anomalous dimension in Eq. (87) we read off 
that the natural scales to evaluate F T are /i ~ |rj_ | ~ 
Aqcd and v ~ p~ ~ Q. By evaluating F T at this scale 
and running it to some common scale (/i, v) for all the 
functions in the factorization theorem, the large loga- 
rithms are summed. 





FIG. 15. One-loop soft function diagrams. Only the Wilson 
lines with gluons attached are shown. 



C. Soft Function 

We now calculate the anomalous dimension of the soft 
function S TT in momentum space. We remind the reader 
that the soft function S 11 = 1 and so does not receive any 
QCD corrections. 

The one-loop contribution from gluon exchange be- 
tween Wilson lines at the same position, e.g. S^(0) and 
S n (0), is the left graph shown in Fig. 15 and is given by 



/ sl = -2i 5 V 2 )(rj( 



47T 

|2fc 3 | 



d d k 
(2n) d 



{-k~ + i0)(fc 2 



7T 



+0(r? u e u ) 



M 2 + i0)(-fc++i0) 



2e 2 2e 



1 



(90) 



where we left out the color factor. We calculated this loop 
integral by first performing the k° integral by contours, 
followed by the k 3 integral and the standard integral. 
We have checked that the result is independent of the iO- 
prescription of the soft Wilson lines (the difference is a 
scaleless integral), implying that the result for this graph 
is the same for incoming quarks or outgoing anti-quarks, 
and vice versa. 

The exchange between Wilson lines at different posi- 
tions, e.g. S^(zi) and S n (0), is the right graph shown in 
Fig. [15] and is given by 

fa = -V(^)V>/^> (k± - rJ (w) 

|2fc 3 |-" 
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1 1 1 
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1 1 
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7T/i 2 (r 2 L / A t 2 ) + M 2 



iO) 



1 

rj 

+O(v°,e ) 



The calculation is very similar, except that this time the 
kx integral is performed using the delta function. In the 



last step we take the limit M 2 — » to isolate the IR 
divergences. 

The RG equations for the soft function are given by 
dS ^T = /d 2 k^(^)S(r ± -k ± ), 

For S TT the color factors for graphs SI and £2 are 

" 4tr[T A T c [T c , T B }} tr [T A T B ] = 2C A , (93) 
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(94) 



Note that the IR divergences again cancel in the sum 
of the diagrams. From Eq. (94) we conclude that the 
natural scales for evaluating the soft function are v ~ 
/i ^ |r_L | ~ Aqcd- An important cross check is provided 

by 



27^ + 7^ 5(l-x 1 )S(l-x 2 )=0, 



(95) 



since the rapidity divergences should cancel between the 
collinear and soft sectors. 



D. Rapidity Resummation 

We now turn to solving the ^-RGE. For definiteness we 
discuss this for the color-correlated soft function, whose 



RGE and anomalous dimension are given in Eqs. (92) 



and ([94]). This is most easily solved in Fourier space, 
where the convolution becomes a product. The anoma- 
lous dimension in position space is given by 



2a s C A 



d 2 r ± e" 



V 2 (r 2 Jv 2 ) + 



2a s C A ,_ z^/Te 



2 P 2 7£ 



In 



(96) 
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where there is no l/(27r) 2 in the Fourier transform since 
there is none in Eq. ( [92] ) . It is now straightforward to 
solve the RGE, since it is local in zj_, 



"^ = ^(.1)^1), 



(97) 



yielding the evolution kernel 



2 /i 2 e 2 ^ 



uj 



2a s (fi)C A ln v 



(98) 



7T 

Transforming back to momentum space, we find 

S TT (r ± , v) = J d 2 k ± uf T (k ± , i/, ^o)^ TT (rx - k ± , ^o) , 
^ e -2 7 ^ r(l-c^) 1 

7T r(i + cj)^(k 2 L ) 1 -^' 1 ; 



^ (kj_,I/,I/ ) = 



This solution is only valid for uj < |, though it can be 
analytically continued for uj < 1. The singularity for 
co> — >• 1 has been observed before in transverse momen- 
tum resummation for Drell-Yan like processes [37] . The 
problem is that Eq. (98) develops a singularity at z± = 0. 



This UV region should not contribute and in Ref. [48] it 
was essentially cut off. We resolve this issue through our 
choice of renormalization scale: Rather than ~ |rj_|, 
we choose vq = 2e~ lE /|zj_|. This scale choice leads to 



Z_l,Z/, 



2e~ lE - 



exp 



a s (/i)C^ 



(L 2 + 2Lln-) 



L = In 



z 2 /i 2 e 2 ^ 



2e- 7s 



Its Fourier transform is given by 

S TT (r ± ,v) = Jdk ± uf T (k ±) 

xS™Yrx-k X) ^) 
V zx / 



1*1 

2e -7E - 



\1, 

(100) 

) (ioi) 



C/r(kx,,,-^) = ^/d Z iJ (|zx||kx|) 



x exp 



7T V \1 



The argument 2e _7s /|zx| is kept in momentum space 
evolution kernel and soft function as a reminder of our 
original scale choice for z/q, though z± is of course no 
longer a variable, since it has been integrated over. In 
the integral for in Eq. ( 101 ), only the region where 

|zj_| ~ 1//^ ~ contributes. We emphasize that the 

scale choice for v§ cancels between the evolution and the 
soft function up to the order in a s that one is working 
at. Our choice for vq rearranges the resummed pertur- 
bation theory, which is required to make the evolution in 
Eq. ([99]) well behaved. 



Our solution is similar to the one proposed in Ref. [49] . 
There an analytic expression for the leading (1 — uj) piece 
of U*j in Eq. (Toil was also obtained, using Borel re- 
summation. 



E. Interference Double PDF 



In Sec. IV Bl we calculated the anomalous dimension for 
the regular dPDF, which we repeat here for the interfer- 
ence dPDF. The results are shown in table |IIl[ Graph // 
contains a matrix in spin space, 



X 



(l 



-1 o\ 
1 
2/ 



(102) 



where the rows and columns correspond to the spin struc- 
tures q, Aq, 5q, respectively. 

Summing up these diagrams, we find that the anoma- 
lous dimension is given by 



2 In— 6(1 

p- 




x) 



2x 



(1-x), 



~-{l-x)X 




(103) 



as a matrix in color space { 1 <g> 1 , T A ®T A } . Repeating the 
analysis around Eq. ( [l4| ) and noting that for N c > 2 the 
matrix multiplying ln(z//p^) has positive eigenvalues, we 
infer that all the interference contributions are Sudakov 
suppressed. 
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I A 



Graph 



Singlet 
0x101 
_ t a t a 



Octet 
\(C A -C d )T A ®T A 

-Cil0l 



7m 



(1-X) + V / p 



lniW5 2 (rx) 



x 1 <g> 1 
T A ®T A 



\{C A + C d )T A ®T A 



Cil ® 1 



[- 



1 _1 1 

7T M 2 ( r 2_/ / x2) + 

-m^^(r ± )] 



75 



Cf 



Cf 



An 



^ + 1)6(1 -x) 



ln^* 2 (r ± ) 



IB' 



-Cf 



-Cf + §CU 



+ ln^^(r ± )] 



7/ 




0x101 



|(CA-C d )T A 0T^ 



-Cil0l 



|(1 -x)X 



W 



C F 



C f 



16(1 -x) 



TABLE III. Diagrams for the interference dPDF. The columns show the graph, the color factors for the 1 1 and T A T A 
color structures (which mix) , the \i and v anomalous dimension and the weight factor. The only spin dependence is in the 
matrix X, given in Eq. (102). 



Interference Soft Function 



In this section we calculate the renormalization of the 
interference soft functions. The graphs are the same as in 
Sec. |VC[ but the color structures are different. We find 
that the the anomalous dimensions of the interference 



soft function are 



7* J (r±) 



Mm) ln i_ § (2) , \ ( SC F -2C f 
tt v 2 y ±J \2C F 10C F -4CV 



6C F 



-2C F 



n 2 n 2 (r 2 ± /n 2 ) + \2C F 10C F - AC A 



(104) 

We will now verify the correctness of this result by 
checking that the interference contribution to the cross 
section is v independent. Ignoring irrelevant factors, 
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C A [4C F S n + (4C A - 12C F )S TT ] 2^[2C F S n + (10C F - AC A )S 




using the anomalous dimensions for the i nterf erence 
dPDFs and soft functions in Eqs. jlQ3| > and ( [105] ) . This 
provides a cross-check of our results. 



VI. CONCLUSIONS 

In this paper we studied double parton scattering from 
the point of view of QCD factorization and renormaliza- 
tion. We presented a detailed derivation of the cross sec- 
tion for double Drell-Yan production in terms of double 
PDFs. Spin and color correlations, as well as interference 
effects, lead to a large number of different contributions. 
In the color-correlated and interference terms, the effects 
of soft radiation are nontrivial, and are given by soft func- 
tions in the factorization formula. We also derived the 
QCD evolution of the dPDFs and soft functions, treat- 
ing both the ultraviolet and rapidity divergences. The 
solution of these equations shows that color-correlated 



and interference contributions are Sudakov suppressed, 
and thus small for double parton scattering at high en- 
ergies. We also discussed several basic properties of dou- 
ble PDFs, such as their classification, properties under 
discrete symmetries and interpretation. In a forthcom- 
ing publication [31], we will discuss the mixing of single 
PDFs into double PDFs and the related issue of double 
counting between single and double parton scattering. 
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